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Abstract 

The Gaussian width of Fermi momentum, p F , is the most important parameter 
of the ACCMM model, and its value is essential in the determination of IV^/V^I 
because the experimental analysis is allowed only at the end-point region of inclu- 
sive semileptonic -B-decay spectrum. We extract the value of \V u b/V c b\ as a function 
of p F . We also calculate the parameter p F in the relativistic quark model using 
the variational method, and obtain p F = 0.54 GeV which is much larger than the 
commonly used value, ~ 0.3 GeV, in experimental analyses. When we use p F = 0.5 
GeV instead of 0.3 GeV, the value of \V u b/V c b\ from ACCMM model is increased 
by a factor 1.81, and can give a good agreement with Isgur et al. model. 
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1. Introduction 



In the minimal standard model CP violation is possible through the CKM mix- 
ing matrix of three families, and it is important to know whether the element V u b 
is non-zero or not accurately. Its knowledge is also necessary to check whether the 
unitarity triangle is closed or not |lj . However, its experimental value is very poorly 
known presently and its better experimental information is urgently required. At 
present, the only experimental method to measure V u b is through the end-point 
lepton energy spectrum of the inclusive 5-meson semileptonic decays, e.g. CLEO 
[|J and ARGUS ||, and their data indicate that V u b is non-zero. Recently it has 
also been suggested that the measurements of hadronic invariant mass spectrum 
H as well as hadronic energy spectrum || in the inclusive B — > X c ^lu decays 
can be useful in extracting \V u b\ with better theoretical understandings. In future 
asymmetric B factories with vertex detector, they will offer alternative ways to 
select b — > u transitions that are much more efficient than selecting the upper end 
of the lepton energy spectrum. 

The simplest model for the semileptonic 5-decay is the spectator model which 
considers the decaying 6-quark in the 5-meson as a free particle. The spectator 
model is usually used with the inclusion of perturbative QCD radiative corrections 
[|J. Then the decay width of the process B — > Xglu is given by 

T B {B - X q lv) ~ T b (b - qlv) = |H,| 2 (^f - \- 9 C^)l (1) 

1927T' 3 nib 6 Tt nib 

where m q is the mass of the final g-quark decayed from 6-quark. As can be seen, 
the decay width of the spectator model depends on mf, therefore small difference 
of nib would change the decay width significantly. 

Altarelli et al. proposed for the inclusive 5-meson semileptonic decays their 
ACCMM model, which incorporates the bound state effect by treating the 6-quark 
as a virtual state particle, thus giving momentum dependence to the 6-quark mass. 
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The virtual state 6-quark mass W is given by 



W 2 (p) =m 2 B + ml 



sp 



2m B \Jv 2 + m 2 s . 



(2) 



in the 5-meson rest frame, where m sp is the spectator quark mass, is the 
B-meson mass, and p is the momentum of the 6-quark inside S-meson. 

For the momentum distribution of the virtual 6-quark, Altarelli et. al. consid- 
ered the Fermi motion inside the 5-meson with the Gaussian momentum distri- 
bution 



where the Gaussian width, p F , is treated as a free parameter. Then the lepton 
energy spectrum of the 5-meson decay is given by 



where p m ax is the maximum kinematically allowed value of p = |p|. The ACCMM 
model, therefore, introduces a new parameter p F for the Gaussian momentum 
distribution of the 6-quark inside 5-meson instead of the 6-quark mass of the 
spectator model. In this way the ACCMM model incorporates the bound state 
effects and reduces the strong dependence on 6-quark mass in the decay width of 
the spectator model. 

The Fermi momentum p F is the most essential parameter of the ACCMM model 
as we see in the above. However, the experimental determination of its value from 
the lepton energy spectrum has been very ambiguous, because various parameters 
of ACCMM model, such as p F , m q and fitted all together from the limited 

region of end-point lepton energy spectrum, and because the perturbative QCD 
corrections are very sensitive in the end-point region of the spectrum. Recently, 
ARGUS || extracted the lepton energy spectrum of B — > X c lu for the whole region 
of electron energy, but with much larger uncertainties. We argue that the value 
p F ~ 0.3 GeV, which has been commonly used in experimental analyses, has no 




(3) 
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theoretical or experimental clear justification, even though there has been recently 
an assertion that the prediction of heavy quark effective theory approach ||, far 
from the end-point region, gives approximately equal shape to the ACCMM model 
with p F ~ 0.3 GeV. Therefore, it is stongly recommended to determine the value 
of p F more reliably and independently, when we think of the importance of its role 
in experimental analyses. It is particularly important in the determination of the 
value of |V^fe/V^fe|, as we explain in section 2. A better determination of p F is also 
interesting theoretically since it has its own physical correspondence related to the 
Fermi motion inside .B-meson. In this context we calculate theoretically the value 
of p F in the relativistic quark model using quantum mechanical variational method 
in section 3. And we obtain p F = 0.54 GeV which is much larger than 0.3 GeV. 
Section 4 contains the conclusion. 



2. Dependence of |V^&/Vd>| on the Fermi momentum parameter p F 



The ACCMM model provides an inclusive lepton energy spectrum of the B- 
meson semileptonic decay to obtain the value of |Ki&/Kf>|- The leptonic energy 
spectrum is useful in separating b — > u transitions from b — > c, since the end-point 
region of the spectrum is completely composed of b — > u decays. In applying this 
method one integrates (01) in the range 2.3 GeV < E\ < 2.6 GeV at the I?-meson 



rest frame, where only b — > u transitions exist [ 10[ . So we theoretically calculate 



/■2.6 dT B 

T{p F )= dEi——(p F ,m sp ,m q ,m B ). (5) 
J2.3 ahji 

In (|) we specified only p F dependence explicitly in the left-hand side, and dYs/dEi 
in the right-hand side is from ([|). Then one compares the theoretically calculated 
r(p F ) with the experimentally measured width T exp in the region 2.3 GeV < E\ < 
2.6 GeV, to extract the value of \V u b\ from the relation 

r exp = \v ub \ 2 x f{p F ) . (6) 



In the real experimental situation 0, ^ |8], [I(J, the only measured quantity is the 
number of events in this region of high E\ compared to the total semileptonic 
events number, i.e. the branching- fraction T exp /Y t ° t l al . Since the value T*°* a ' is 
propotional to |Kb| 2 , only the combination \V ub /V cb \ 2 is extracted. 

We now consider the possible dependence of \V u b/V c b\ 2 as a function of the 
parameter p F from the following relation 
f 



exp 

T^total 
1 s.l. 



V ub 2 



V, 



do 



xT( P[ 



p F =p F 



V ub 2 



Vr 



cb 



x T(p F = 0.3) 



(7) 



Pf 



=0.3 



where \V ub /Vcb\p F = PF is determined with an arbitrary value of the Fermi momentum 
parameter p F . In the right-hand side we used p F =0.3 GeV because this value is 
commonly used in the experimental determination of IK^/V^. Then one can get 
a relation 
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v ub 


v cb 


P F =P F 


v cb 



p„=0.3 



£(03) 

Wf) 



(8) 



In section 3, we obtain p F = 0.54 GeV using the variational method in the 
relativistic quark model. If we use p F = 0.5 GeV, instead of p F = 0.3 GeV, in the 
experimental analysis of the end-point region of lepton energy spectrum, the value 
of |V^6/V^6| becomes significantly changed. We numerically calculated theoretical 



ratio r(0.3)/r(0.5) by using @ and (§ with 



m 



0.15 GeV, m q = 0.15 GeV, 



which are the values commonly used by experimentalists, and = 5.28 GeV, to 
get its value as 1.81, which finally gives 
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v ub 


v cb 


p F =0.5 


v cb 



x 1.81 . 



(9) 



p„=0.3 



Previously the CLEO [HJ analyzed with p F = 0.3 GeV the end-point lepton 
energy spectrum to get 

10 2 x \V ub /V cb \ 2 = 0.57 ±0.11 (ACCMM [7]) 

= 1.02 ±0.20 (Isgur et.al. 0) . (10) 

As can be seen, those values are in large disagreement. However, if we use p F = 0.5 
GeV, the result of the ACCMM model becomes 1.03, and these two models are 



in a good agreement for the value of |Vu&/Vd,|. Finally we show the values of 
\Vub{p F )/V ub (p F = 0.3)| as a function of p F in Fig. 1. 

3. Calculation of p F in the relativistic quark model 

We consider the Gaussian probability distribution function 0(p) in (3) as the 
absolute square of the momentum space wave function x(p) °f the bound state 
.B-meson, that is, 

0(p) = 47T| X (p)| a , X(P) = ( ^ )3/2 e- p2/2 ^- (11) 

The Fourier transform of x(p) gives the coordinate space wave function ip(r), which 
is also Gaussian, 

M T ) = (^)3/2 e -rV F /2_ (12) 
'7T 



Then we can approach the determination of p F in the framework of quantum 
mechanics. For the 5-meson system we treat the 6-quark non-relativistically, but 
the u- or <i-quark relativistically with the Hamiltonian 



H = M + M + V / P 2 + m2 + ( 13 ) 
where M = ra& is the 6-quark mass and m = m sp is the u- or <i-quark mass. We 
apply the variational method to the Hamiltonian ( |i3|) with the trial wave function 

^(r) = (^f/V^ 2 / 2 , (14) 



where \x is the variational parameter. The ground state is given by minimizing the 
expectation value of H, 

(H) = = E(fi), ^E(fi) = at fM = jl, (15) 

and then p, = p F and E = E{p) approximates tub- The value of \i or p F corre- 
sponds to the measure of the radius of the two body bound state as can be seen 
from (r) = -?=- and (r 2 )^ = |-. 



In (|T^) we take the Cornell potential which is composed of the Coulomb and 
linear potentials, 

rv 

V(r) = -— + Kr. (16) 
r 

For the values of the parameters a c (= |a s ), K, and the 6-quark mass rrib, we use 
the values given by Hagiwara et al. fll2|, 



a c = 0.47 (a s = 0.35), K = 0.19 GeV 2 , m b = 4.75 GeV, (17) 

which have been determined by the best fit of all the known (cc) and (bb) bound 
states. For comparison we will also consider a c = 0.32 (a s = 0.24), which corre- 
sponds to a s (Q 2 = m%). 

Before applying our variational method with the Gaussian trial wave function 
to the S-meson system, let us check the method by considering the T(bb) system. 
The Hamiltonian of the T(bb) system can be approximated by the non-relativistic 
Hamiltonian 

P 2 

H~2m b + — + V(r). (18) 
m b 

With the parameters in (|T7]) (or with a c = 0.32), our variational method with 
the Gaussian trial wave function (14) gives p F = p, = 1.1 GeV and E = E(p) = 



9.49 GeV. Here p F = 1.1 GeV corresponds to the radius R(T) = 0.2 fm, and 
E(T) = 9.49 GeV is within 0.3 % error compared with the experimental value 
E exp = my = 9.46 GeV. Therefore, the variational method with the non-relativistic 
Hamiltonian fll8|) gives fairly accurate results for the T ground state. 

However, since the u- or d-quark in the 5-meson is very light, the non-relativistic 
description can not be applied to the S-meson system. For example, when we ap- 
ply the variational method with the non-relativistic Hamiltonian to the S-meson, 
we get the results 

p p = 0.29 GeV, E = 5.92 GeV for a s = 0.35, (19) 
p F = 0.29 GeV, E = 5.97 GeV for a s = 0.24. (20) 



7 



The above masses E are much larger compared to the experimental value m# = 
5.28 GeV, and moreover the expectation values of the higher terms in the non- 
relativistic perturbative expansion are bigger than those of the lower terms. There- 
fore, we can not apply the variational method with the non-relativistic Hamiltonian 
to the 5-meson system. 



Let us come back to our Hamiltonian (|13|) of the .B-meson system. In our varia- 
tional method the trial wave function is Gaussian both in the coordinate space and 
in the momentum space, so the expectation value of H can be calculated in either 
space from (H) = (il>(r)\H\ip(r)) = (x(p)\H\x(p)) ■ Also, the Gaussian function 
is a smooth function and its derivative of any order is square integrable, thus any 
power of the Laplacian operator V 2 is a hermitian operator at least under Gaussian 



functions. Therefore, analyzing the Hamiltonian ([13]) with the variational method 
can be considered as reasonable even though solving the eigenvalue equation of 
the differential operator flT3| ) may be confronted with the mathematical difficulties 
because of the square root operator in (JT3|) . 

With the Gaussian trial wave function ( P"T|) or (13), the expectation value of 
the Hamiltonian (|i~3|) can be calculated easily besides the square root operator, 



(p 2 ) = Mr)|p 2 |^(r)) = ( X (p)|p 2 | X (p)) = ^ 2 , (21) 
(V(r)) = (^r)\-^ + Kr\^(r)) = —(-a c v + K/ri. (22) 



Now let us consider the expectation value of the square root operator in the mo- 
mentum space 

p 2 + m 2 > = (x(p)| V / P 2 + ^ 2 |x(p)) = (4=) 3 J™ e~ p2/fl2 VP 2 + m 2 d 3 p 
= % e- x \lx 2 + (m//x) 2 x 2 dx. (23) 

V7T JO V 



The integral ( PB| ) can be given as a series expansion by the following procedure. 
First, define 

roo poo „ 

I(s) = / Vx^Tli x 2 e~ x dx = s 2 / V¥Tlt 2 e- st dt, (24) 
Jo Jo 



I (s) = / V^TTse- x dx = s VWTl e~ st dt. 



(25) 



Next, from (|24] ) and (p5|), we find the following differential relations 



Combining two equations in (p6|), we get a second order differential equation for 

sl"(s) - (i + s)I'(s) + 1/( S ) = 0. (27) 
The series solution to fl27|) is given as 



1(a) 



ciJi(s) + c 2 / 2 (s) 



F(s; -, 3) = s 2 { 1 + -s + —s 2 + — s 3 + —s A + ■■•}, (28) 
v ' 2' ; I 2 32 192 1024 J' v ; 



h{s) 



- ^1 



i T^rds = s 2 In s( 1 H — s -\ s 2 

[his)} 2 16 V 2 32 

5 o To 7 a 

s H s 2 H s 3 H s 4 + ■ 

2 32 192 1536 



where F(s; |,3) is the confluent hypergeometric function which is convergent for 
any finite s, and the integral constants c\ ~ —0.095, c 2 = — 1. See Appendix for 
the derivation of these numerical values for q. 

Finally, collecting (|2lf) , (|22[) and (|2"3f) , the expectation value of if is written as 



2fl 



1 + hm/tf + (A - 2 Cl )(m//i) 4 + i(m/^) 4 ln(m//i) 



, (29) 



up to (m/nY 



With the input value of m = m sp = 0.15 GeV, we minimize (H) of (|29"D, and 
then we obtain 



Pf = /i = 0.54 GeV, ffig = E = 5.54 GeV for a s = 0.35, (30) 
ft = 0.49 GeV, E = 5.63 GeV for a s = 0.24. 



9 



Here let us check how much sensitive our calculation of p F is by considering the 
case where m = m sp = for comparison. For m sp = the integral in (^) is 
done easily and we obtain the following values of // = p F by the above variational 
method. 

fl = 0.53 GeV, E = 5.52 GeV for a s = 0.35, (31) 
ft = 0.48 GeV, E = 5.60 GeV for a s = 0.24. 

As we see in (j3~i~D , the results are similar to those in ( j30|) where m sp = 0.15 GeV. 
We could expect this insensitivity of the value of p F to that of m sp because the 
value of m sp , which should be small in any case, can not affect the integral in (23) 
significantly. 

The calculated values of the 5-meson mass, E, are much larger than the mea- 
sured value of 5.28 GeV. The large values for the mass are originated partly because 
the Hamiltonian ( |29| ) does not take care of the correct spin dependences for B and 
B* . The difference between the pseudoscalar meson and the vector meson is given 
arise to by the chromomagnetic hyperfine splitting, which is given by 

K = ^5W 2 V 2 (--)- (32) 
SMm r 

Then the expectation values of V s are given by 

and we treat (V^) only as a perturbation. Then with the input value of m = m sp = 
0.15 GeV, we get for B meson 

p F = 0.54 GeV, E B = 5.42 GeV for a s = 0.35, (34) 
p F = 0.49 GeV, E B = 5.56 GeV for a s = 0.24, 

p F = 0.54 GeV, E B * = 5.58 GeV for a s = 0.35, (35) 
p F = 0.49 GeV, E B * = 5.65 GeV for a s = 0.24. 
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and for B* 



The calculated values of the £?-meson mass, 5.42 GeV (a s = 0.35) and 5.56 
GeV (a s = 0.24) are in reasonable agreement compared to the experimental value 
of tub = 5.28 GeV; the relative errors are 2.7% and 5.3%, respectively. However, 
for the Fermi momentum p F , the calculated values, 0.54 GeV (a s = 0.35) and 0.49 
GeV (a s = 0.24), are larger than the value 0.3 GeV, which has been commonly 
used in the experimental analyses of energy spectrum of semileptonic _B-meson 
decay. The value p F = 0.3 GeV corresponds to the B-meson radius Rb ~ 0.66 fm, 
which seems too large. On the other hand, the value p F = 0.5 GeV corresponds 
to Rb ~ 0.39 fm, which looks in reasonable range. 

4. Conclusion 

The Gaussian width of Fermi motion, p F , is the most important parameter of 
the ACCMM model, and the value p F ~ 0.3 GeV has been commonly used in ex- 
perimental analyses without clear theoretical or experimental evidence. Therefore, 
it is recommended to determine the value of p F more reliably, when we think of 
its importance in experimental analyses. We calculated the value for p F in the rel- 
ativists quark model using the variational method. We obtained p F = 0.54 GeV, 
which is much larger than 0.3 GeV. We also derived the ground state eigenvalue 
of E B — 5.5 GeV, which is in reasonable agreement with the experimental value 
m B = 5.28 GeV. 

We studied the dependence of \V u b/V c b\ on the Fermi momentum parameter 
p F in the ACCMM model, and extracted \V ub /V cb \ as a function of p F . It shows 
that |T4b/Kb| is very much dependent on the value of p F . When we use p F = 0.5 
GeV instead of 0.3 GeV, (Kb/V^l is increased by a factor 1.81. Then the previous 
discrepancy between the ACCMM model and the Isgur et al. model for the value 
of 1 14b /K& | 2 turns into a good agreement. 
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Appendix 

The integration constants c\ and c 2 in (p8|) are given by the following relations, 



7(0) 



1 f°° 3 _ x 2 1 

-c 2 = / are da; = -. 

2 Jo 2 



2"(a«0) = 2c 1 + c 2 (-^ln S -^; 



1 Z" 00 2 

— / x 2 (x 2 + sy 3/2 e~ x dx at s « 0. 
4 Jo 



Then, from (|36|) , we get 



C2 



-1. 



The integral in (07|) can be expanded as 

J(s = a 2 ) = r x 2 (x 2 + a 2 )-^ 2 e- x2 dx 
Jo 

= I" x 2 [(x + a) 2 -2axr 3/2 e- x2 dx 
Jo 





oo 

E 



x 2 (x + a) 3 



(2n + l)!a r ' 



2ax 



-i -3/2 



(x + a) 2 



e x dx 



n =o 2«(n!) 2 Jo (x + a) 2n + 3 
Next the integral in (|39"D is obtained by 



e x dx. 



/ 

Jo 



X n+2 „2 . 1 



-e x dx 



f) \ 2n ~*~ 2 roc T n+2 

£ / e~ x dx. (40) 

oa / Jo x + a 



(x + a) 2ri + 3 ' (2n + 2)! \da J Jo .m « 
Again the integral in ( EH ) is related to another integral 



.oo x n+2 ^ n+1 (_ (/ ,/V ( _ ; , > 

Jo x + a 



\n+2 



(36) 



(37) 



(38 



(39) 



/"GO p x 

/ dx. (41) 

Jo x + a 



The integral in ( f4l|) is given, for an infinitesimal value of a, by integration by parts, 

_ 2 f°° _ 2 7 

dx = e x /n(x + a)| XD — / e x (— 2x)ln(x + a)dx = — In a 1- 0(a), 

Jo 2 

(42) 



o x + a 



where 7 ~ 0.5772 is the Euler's constant. Collecting fl40|), ©, and (||), 

1 / ft \ 2n+2 



7 



— p-ry ra n (^l (-a) rt+2 {-/n a- - + 0(a)}. (43) 

^2"(n!) 2 (2n + 2) / 2 
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To get the constant c±, we should extract a logarithmic term and constants from 

o oo /_-. Nn+l / Q \ 2n+2 

v ; v 2 2 ; ^ 2»(n!) 2 (2n + 2) V<9a/ v ; ' v ' 



n n ( d Y n+2 (n n + 2 ]r] n){ "f 2 (2^ + 2)! (n + 2)!(2n + l-fc)!(-l) fc + 1 
a ^-j (a lna)\ a=0 - ^ ^- - - - (n + 2 _ fc)! 

= (-l) n+1 (n + 2)!(n- 1)!. (45) 



Inserting (|45|) into (0), we get 



J(a « 0) = -Zn a - | - 1 + /?, (46) 



where (3 = ££° =1 1/(^2") « 0.6932. Then, from (|7|), we get 



^^-s^-^"- - 0975 ' (47) 
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Fig.l The ratio \V u b(p F )/V u b(p F = 0.3) | as a function of p F . 
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